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ABSTRACT 

Variable stars have been identified for tlie first time in tlie very metal-poor Blue Compact dwarf 
galaxy IZwl8, using deep multi-band (F606W, F814W) time-series photometry obtained with the 
Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). We detected 34 
candidate variable stars in the galaxy. We classify three of them as Classical Cepheids, with periods 
of 8.71, 125.0 and 130.3 days, respectively, and other two as long period variables with periodicities 
longer than a hundred days. These are the lowest metallicity Classical Cepheids known so far, thus 
providing the opportunity to explore and fit models of stellar pulsation for Classical Cepheids at 
previously inaccessible metallicities. The period distribution of the confirmed Cepheids is markedly 
different from what is seen in other nearby galaxies, which is likely related to the star bursting nature 
of IZwl8. The long period Cepheids we have detected in IZwl8 seem to indicate that massive stars 
at the metallicity of IZwl8 (Z=0.0004) may cross the instability strip long enough to be observed. 

By applying to the 8.71 days Cepheid theoretical Wesenheit (V,/) relations based on new pulsation 
models of Classical Cepheids specifically computed for the extremely low metallicity of this galaxy 
(Z=0.0004, Y=0.24), we estimate the distance modulus of IZwl8 to be ^XQ= 31.4 ± 0.3 (D=19.0t2:5 
Mpc) for canonical models of Classical Cepheids, and of 31.2 ± 0.3 mag (0=17.4^2 2 Mpc) using over 
luminous models. The theoretical modeling of the star's light curves provides /xq = 31. 4± 0.2 mag, 
D=19.0l|^'7 Mpc, in good agreement with the results from the theoretical Wesenheit relations. These 
pulsation distances bracket the distance of 18.2 ±1.5 Mpc inferred bv lAloisi et al.l (|2007l ) using the 
galaxy's Red Giant Branch Tip. 

Subject headings: — galaxies: individual (IZwl8) — stars: variables: other - Classical Cepheids — 
distance scale — galaxies: Blue Dwarf Galaxy — techniques: photometry 



1. INTRODUCTION 

Discovered by iSearle et al.l ()1973f ) more than thirty 
years ago the Blue Compact Dwarf galaxy IZwl8 
remains one of the mos t metal poor (1/50 Zq, 
iSkillman fc Kennicutti 119930 and intriguing objects in 
the local Universe. The galaxy is very rich in gas content 
and blue stars. These two features both suggested that 
IZwl8 might be undergoing its first episode of star for- 
mation, thus representing one of the best local analogs 
to primordial galaxies in the distant Universe. 

Several HST studies by different authors have tried 
to characterize the evoluti onary status of IZ wl8, with 
rather discordant results. lAloisi et al.l ([1999') were the 
first to detect Asymptotic Giant Branch (AGB) stars in 
the galaxy, and thus to demonstrate that IZwl8 is as old 
as several hundreds Myr, at least, using HST Wide Field 
Plane tary Camera 2 (WFPC2) data. Ilzotov fc ThuanI 
()2004[ ) failed to detect red giant branch (RGB) stars in 
IZwl8, based on HST/ACS observa tions of the galaxy. 
However, PMomanv et al.l (j2005D and iTosil (|2006( ). based 

* BASED ON OBSERVATIONS WITH THE NASA/ESA HUB- 
BLE SPACE TELESCOPE, OBTAINED AT THE SPACE TELE- 
SCOPE SCIENCE INSTITUTE, WHICH IS OPERATED BY 
AURA INC., UNDER NASA CONTRACT NAS 5-26555. 



on a re-analysis of the same ACS data set, suggested 
that IZwl8 is older than previously believed, since they 
found some evidence for the galaxy's RGB Tip around 
/ - 27 mag. Our new time-series HST/ACS data (HST 
program Id. 10586, PI: Aloisi) have allowed us to defini- 
tively solve the mystery on the age of IZwl8, by pinning 
down the galaxy's distance, using Classical Cepheids that 
we detected for the first time in this galaxy. The dis- 
tance to IZwl8 inferred from the Classical Cepheids is 
in very good agreement with the independent estimate 
obtained from the galaxy's red giant branch (RGB) Tip 
D = 18.2± 1.5 Mpc), which is identified bv lAloisi et"~al] 
2001 at Jo — 27.27 ±0.14 mag (for an assumed Galactic 
foreground extinction of _E(_B— V^)=0.032 mag, and shows 
that IZwl8 is further away than previously thought. In 
this paper we present in detail the variable stars we have 
identified in IZwl8 and discuss in particular the galaxy's 
Classical Cepheids. Observations and data reduction are 
described in Section 2. The techniques applied to iden- 
tify candidate variable stars and derive their periodicities 
are described in Section 3. The properties of the IZwl8 
confirmed variable stars are presented in Section 4, while 
the distance inferred from the Classical Cepheids is pre- 
sented in Section 5. Finally, our results and conclusions 
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are summarized in Section 6. 

2. OBSERVATIONS AND DATA REDUCTIONS 

Time series imaging of IZwl8 in the F606W (broad 
V) and F814W (/) filters was obtained with the Wide 
Field Channel of the Advanced Camera for Surveys 
(ACS/WFC) on board the HST in 13 different epochs 
properly spread over a time interval of 96 days from 2005 
October to 2006 January. The time window and the 
scheduling of the single epoch observations were specifi- 
cally devised to allow the discovery and characterization 
of Classical Cepheids with periods shorter than 100 days, 
if present. 

Each epoch consisted of four optimally dithered single 
exposures per filter, that we corrected for geometric dis- 
tortion and co-added by using the MultiDrizzle software, 
to produce single images per epoch per filter, totaling 
about 35-37 minute integration time per epoch, and to- 
tal integration times, over the 13 epochs, of ~ 7.7 and 
~ 7.3 hours in F606W and F814W, respectively. The 
proprietary data were complemented by archival F555W 
{V) and F814W ACS/WFC images of IZwl8 from HST 
program Id. 9400 (PI Thuan), obtained in 2003, May- 
June, over a time span of 11 days, thus gaining three and 
five additional epochs in the F814W and F555W filters, 
respectively. The additional data allowed us to increase 
the number of available phase-points on the light curves 
to 15 in F814W, 13 in F606W and five in F555W and 
to extend the time interval beyond the 96 days spanned 
by the proprietary observations. This is valuable to pin 
down periods longer than the time window covered by 
the proprietary data, and to reveal long period variabil- 
ity. The log of the observations for both proprietary and 
archival data is provided in Table 1 . 

The severe crowding of the IZwl8 images, the demand 
of achieving the maximum possible depth in the color 
magnitude diagram (CMD), and the need of processing 
the single epoch data to generate light curves for the 
candidate variables, made the photometric reduction of 
our data very challenging and time-consuming, and the 
use of PSF-fitting photometry mandatory. To produce 
the galaxy CMD, we built four coadded images by stack- 
ing all single epoch images per filter for both archival 
and proprietary data, separately. Then we performed 
PSF photometry of each of these four coa dded images us - 
ing the DAOPHOT/ALLSTAR package (Stetson 198^. 
The resuhing CMD was pubhshed by Aloisi ct al. (2007). 
It contains ~ 2100 sources having detection in all four 
deep images {V and / for both datasets), and shows the 
galaxy's RGB Tip at /o==27.27 ± 0.14 mag. 

To produce light curves for the candidate variable ob- 
jects we followed the approach described in the Hq cook- 
book by Turner (1997) and used the MONTAGE package 
within DA OPHOTII/ALLSTAR/ALLFRAME (.Stetson 
119871 119941) to build a ma ster frame out of th e four coad- 
ded images obtained by lAloisi et al.l ()2007[ ). Then we 
performed DAOPHOTII/ALLSTAR/ALLFRAME PSF 
photometry on each single epoch (15 I and 18 V im- 
ages) using as a reference the catalogue of sources ob- 
tained from the master frame. The new photometry 
was calibrated to the F555W, F606W and F814W HST 
Vegamag system b y using the I Zwl8 photometric cata- 
logue presented in lAloisi et al.l ()2007D as a list of stan- 
dard reference stars. Finally, the calibration to the 



Johnson-Cousins photometric system was ob tained using 
the transformations given by iSirianni et all §005) ( see 
lAloisi et al.|[2007l , for details). The F606W and F555W 
photometry were calibrated to the Johnson- Cousins pho- 
tometric system separately, and the resulting standard V 
magnitudes compared. The result of this comparison is 
shown in Figure [TJ for V — I colors redder than 0.4 mag, 
which roughly corresponds to the average color inside the 
classical instability strip (IS) in IZwl8. 

A zero point shift was found to exist between the 
two V^-calibrated photometry for magnitudes brighter 
than V=26.5 mag. (see arrow in Figure In fact, 

the V magnitudes derived from the calibration of the 
F606W data {Vpeoew) are on average 0.04 mag brighter 
than the V magnitudes derived from the calibration 
of the F555W data (Fpsssw)- Since the F555W fil- 
ter more closely resembles the photometric passband 
of the Johnson-Cousins V photometry, and the F606W 
fluxes might be somewhat overestimated because of 
non-negligible gas contamination, the t^reoew magni- 
tudes were corrected for this zero point and tied to 
the Vf555W ones. The photometric reductions with 
DAOPHOTII/ALLSTAR/ALLFRAME produced a new 
I,V ~ I CMD of IZwl8 which contains about 7000 
sources, all having photometric errors ayj < 0.5 mag 
and SHARP parameter^ in the range from —1 to +1, 
and reaches V ~ 30.0 mag, i.e. about half a magnitude 
fainter than the lAloisi etahi (2007) CMD. The new CMD 
is shown in Figure [21 We note that a fraction of the ad- 
ditional stars in our CMD populate the blue plume of 
IZwl8, and that some of them appear to be bluer than 
the Padua tracks. An eye inspection of a map of IZwl8 
confirms that these are likely gas contaminated stars, 
usually located in gas-rich regions and/or gas filaments. 

3. VARIABLE STAR IDENTIFICATION AND PERIOD 
SEARCH 

We used four independent approaches to identify vari- 
able stars in IZwl8 and to select bona fide variables 
from false detections. First we applied to our proprietary 
F606W and F814W data sets the Optical Image Subtrac- 
tion Techniq ue as performed within the package ISIS 2.2 
(|Alardll2000l ). This method has shown to be very pow- 
erful to detect variables in highly crowd ed stellar fields 
such as the main body of IZwl8 (see e.g.. lBaldacci et all 
120051 ) . With this approach we obtained a list of candidate 
variables and differential fiux light curves, on which to 
perform the period search. However, the under-sampling 
of the ACS/WFC images revealed a very challenging is- 
sue for ISIS. The number of sources flagged as candidate 
variables was extraordinary large, making the check of 
the individual objects extremely time-consuming and dif- 
ficult. 

Moreover, inspection of the images revealed that ISIS 
had very often failed to distinguish real stars from spu- 
rious sources, galaxies and background noise. 

We then used the new 

DAOPHOTII / ALLSTAR / ALLFRAME photometry 
of the single-epoch data and an output parameter of 
ALLFRAME, the Welch-Stetson variability index (see 

^ The SHARP parameter is related to the intrinsic angular size 
of the object image and measures the regularity and symmetry of 
the PSF stellar profile. 
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iWelch fc StetsonI 119931 and the cookbook) to flag 
candidate variable stars in the catalogue of measured 
sources returned by ALLFRAME. This second approach 
allowed us to pick up 66 candidate variables. All of them 
were already present in the list of variables detected by 
ISIS. However, by using the light curves in magnitude 
scales produced by ALLFRAME, it was much easier and 
straightforward to check the candidate variables than 
by working with the differential flux light curves of ISIS, 
where we lacked information on the actual amplitude of 
the light variation. 

As a third approach, PSF-fitting photometry of 
the single epoch dat a was performed with DoPHOT 
(jSchechter et al.l 119931 ) and we used the pack age proce- 
dures, as discussed in iSaha fc Hoessell ()199Cl[ ). to search 
for periodic variables. Briefly, once we had star magni- 
tudes at several epochs, variable star candidates were de- 
tected as the sources whose data points had a large scat- 
ter in the A MAG vs MAG plane (where A MAG is the 
difference between magnitudes at two different epochs) 
and good parameters of sharpness and goodness-fit cri- 
teria (x^). To make the flag as candidate variable more 
robust, the empirical scatter was compared with the av- 
erage error estimated by DoPHOT. This procedure was 
applied to several frame pairs. The scattering data points 
were then checked by eye. They often revealed spurious 
variable sources, due to the difficulty of performing a re- 
liable PSF-fitting photometry in IZwl8. 

Finally, identification of candidate sources was also 
performed by visual inspection of the F606W and 
F814W images by one of us (A.S.). In this approach each 
pixel is the x^ value calculated from the corresponding 
pixels in the co-registered target images at all epochs. 
The stability of the HST/ACS PSF, combined with the 
careful registration of the images makes it possible to 
visualize the variable objects very effectively. The co- 
registered images in the same passband can be stacked 
as a data cube, with each layer in the "z" direction as one 
of the time epochs. An estimate of the noise in each pixel 
of this data cube is obtained from the Poisson statistics 
and the read noise of the detectors. From the variance 
of each pixel along the "z" direction, and the noise es- 
timates as above, a x^ is obtained, thus producing a x^ 
image. All of the definite variables identified in this paper 
are clearly visible as pixel groupings with high values in 
the x^ images. We have not attempted a detailed PSF 
analysis of these x^ "stars" , but visually, the majority 
of these enhancements appear as compact and "round" 
structures, with sizes consistent with that of single stars. 
Note that the level of enhancement corresponds to the 
statistical significance for variability, and objects that 
are relatively very faint can show up prominently. Note 
also that several x^ enhancements which are consistent 
with the "footprint" of single stars, appear in the central 
body of IZwlS. In these areas, the crowding makes it 
impossible to perform any photometry, or even de-blend 
individual stars in the actual images. It is thus unfortu- 
nate that even though the x^ maps show the presence of 
variable stars so cleanly in such confused regions, they 
remain inaccessible to further analysis. The x^ image ap- 
proach is an extension of the "difference image" method 
for identifying variables. It would eliminate the need to 
compare images pair-wise, allowing rather to look for sta- 



tistically significant variability in multi-epoch data, while 
preserving the information provided by the correlation 
between neighboring pixels due to the PSF. Clearly, the 
method bears further investigation and development. 

Candidate variable stars found by the four indepen- 
dent procedures described above were counter identified 
by coordinates, returning a final catalogue of 34 bona 
fide candidates showing light variations larger than three 
times the photometric error at their average magnitude, 
that we retained for further analysis. Nineteen out of 
the 34 candidates lie on and around the main body of 
IZwlS, the remaining 15 ones on the galaxy secondary 
body. Figure [3] shows the location of the 34 candidate 
variables on a map of IZwlS. 

Period search and study of the light curves for the bona 
fide candidate variables were performed using an itera- 
tive procedure that we applied to the FS14W, F606W 
and F555W instrumental magnitudes, separately. A first 
guess of the star's periodicity was obtained using the 
Phase Dispersion Minimization (PDM) algorithm, within 
the IRAF environment, making a quick search over a 
large time interval ranging from 2-3 days to about 200 
days. The period refinement was then obtained with the 
Graphical Analyzer of Time Series (GRaTiS), a private 
software develop ed at the Bologna Observatory (see e.g., 
iClementini et al . 2000), which uses b oth the Lomb peri- 
odogram (|Lombil976t lScargleiri9S2[ ) and the b est fit of 
the d ata with a truncated Fourier series (e.g., iBarningI 
[1961 . First we run GRaTis on the 15 epochs of the 
2003-^2007 FS14W data set deriving a period that we 
used to phase the F555W and F606W light curves, and 
iteratively improved by comparing the light curves in the 
three different pass-bands. Candidate variables showing 
the same consistent periodicity in all the three bands 
were confirmed as bona fide variable stars, while the 
other stars remained potential candidates. This proce- 
dure allowed us to select, among the 34 candidates, 5 
variable stars with consistent light curves in all the three 
filters. These five bona fide variables are listed in Table 2 
and plotted as red filled triangles in Figure [2] Their po- 
sition on the body of IZwlS is shown in Figure H) Four of 
them (VI, V4, V6 and V7) lie on the galaxy main body. 
The fifth one, V15, is on a gas shell surrounding the 
main body and pointing towards the direction of IZwlS 
secondary body. The 5 confirmed variables are further 
discussed in Section 4. 

For the remaining 29 objects (see Table 3) we only have 
an identification, the mean V , I magnitudes provided by 
ALLFRAME, and a rough indication of the light vari- 
ation in the F606W and FS14W bands, since several of 
these objects lie in areas where crowding makes it impos- 
sible to perform a reliable photometry of the single epoch 
data. We also note that although many of the objects 
in Table 3 show significant multi band light variations, 
for most of them the V light curve does not match the / 
curve for the same period. This is likely due to crowding 
that makes the counter identification of an object in dif- 
ferent photometric bands rather uncertain, and to stellar 
blending, which affects in different ways the star magni- 
tude, depending on the photometric band. It could also 
be possible that some of these candidate variables are in 
unresolved stellar associations. The candidate variable 
stars are shown as orange filled circles in Figure [2] 
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4. IZW18 CONFIRMED VARIABLE STARS 

The time series data of the five bona fide variable stars 
were cahbrated to the Johnson-Cousins photometric sys- 
tem following the general procedure described in Section 
2 and properly taking into account the color variation of 
the variable stars during the pulsation cycle. In practice, 
we associated to each phased F555W and F606W data 
point the F814W value read, at the same phase, from the 
model best fitting the F8I4W data, and to each phased 
F814W data point the F606W value at same phase read 
from the model best fitting the F606W data. The V 
magnitudes derived from the F606W photometry were 
corrected for the 0.04 mag zero point shift before com- 
bining them to the F555W calibrated photometry to de- 
rive the average magnitudes of the variable stars over 
the pulsation cycle. The time-series V, / photometry of 
the 5 confirmed variables is provided in Table 4 which is 
published in its entirety in the electronic edition of the 
journal. 

Results from the analysis of the 5 confirmed variables 
are summarized in Table 2, where for each star we pro- 
vide identification number, coordinates, period, time of 
maximum light of the V data, intensity-averaged mean 
magnitudes and colors, amplitudes in the V and / bands 
of the Johnson-Cousins photometric system, and ampli- 
tude ratios. Figure [5] shows the Johnson-Cousins V, / 
light curves of the five confirmed variables, where data 
have been folded according to the periods and times of 
maximum light listed in Table 2. Error bars show the 
intrinsic errors of the PSF photometry, as computed by 
ALLFRAME, on the single epoch measurements. For 
the faintest variable star (V6) they are of the order of 
0.09 mag in V and 0.12 mag in /. We note that only 
for V6 the pulsation cycle is repetitively covered by our 
observations, thus making the period determination ro- 
bust. We estimate an uncertainty of the order of ± 0.08 
days on the period of this star. The other 4 confirmed 
variables all have long periodicities exceeding the time 
window covered by the proprietary data. Although the 
addition of the archival data allowed us to improve the 
period determination, for these stars we find more un- 
certain periods than for V6. 

For the longer period variables we first used the Lomb 
algorithm within GRaTiS to explore a period interval 
ranging from a few days up to 300-500 days, and chose, as 
usual, the most likely period as that corresponding to the 
highest peak in the star's power spectrum. The period 
definition was then refined using the truncated Fourier 
algorithm. In general, given the smooth and sinusoidal 
shape of the light curves, 2 harmonics were sufficient to 
best fit the data. According to the width of the highest 
peaks in the power spectra of VI and V15, the accuracy 
of the derived periodicities is of the order of ± Iday. 
However, it should be noted that both these two stars 
also have lower peaks in the power spectra corresponding 
to P=113 and 152 days for VI, and P=110 and 146 days 
for VI 5. We are inclined to attribute these alternative 
periodicities to the windowing of the data. However, only 
further time-series covering much longer time windows 
could allow us to constrain the actual periods. 

The period definition is more uncertain for V4 and 
V7. Indeed, the Lomb periodogram would favor a period 
around 165 days for V7, while the truncated Fourier se- 



ries favors a shorter period of 106 days. Our data do not 
allow us to distinguish between these two periodicities. 
For practical purposes we adopted the shorter period to 
fold the light curve data, but this period should be taken 
with caution. 

Similarly, the very long period of V4 (196 days) signif- 
icantly exceeds the time window covered by the propri- 
etary data (96 days), and after addition of the archival 
data yet we apparently covered only the rising branch of 
the star light curve, thus rendering the period determi- 
nation very uncertain. The highest peak in the power 
spectrum of V4 is double peaked, with two equally pos- 
sible periodicities of 187 and 196 days. Moreover, several 
other longer periods exceeding the 200 days are also pos- 
sible. We used the 196 period to fold the ligh curve data; 
however, as with V7, the period of V4 should be consid- 
ered only a tentative estimate. 

In order to classify the bona fide variable stars in types 
we have used the following three diagnostics: (i) the 
star's pulsation period; (ii) the shape, the morphology, 
and the amplitude ratios of the light curves; (iii) the 
location on the galaxy's CMD in comparison with evolu- 
tionary tracks by the Padua group (Fagotto et al. 199^ 
; hereafter Padua94) for the mass range from 5 to 60 
M© (solid curves in Figure H]). The Cepheid candidates 
were also compared with the theoretical boundaries of 
the classical instability strip (IS) specifically computed 
for the extremely low metallicity of IZwl8 (Z=0.0004, 
Y=0.24, Marconi et al. 2010 accepted on ApJ, hereafter 
MIO), and for two different choices of the mixing length 
parameter, namely, a=1.5 and 2 (solid and dashed lines 
in Figure [5]). 

Three out of the five variable stars, namely VI, V6 
and VI 5, have light curves closely matching those of 
Classical Cepheids. Figure [2] shows the location of the 
five variable stars on the /, (V - I) CMD of IZwl8 ob- 
tained by the present analysis. The confirmed variable 
stars (filled triangles in Figure [2]) are plotted according 
to their intensity-averaged mean magnitudes and colors, 
obtained by weighting over the star pulsation cycle and, 
with the exception of V4 and V7, occupy a region of the 
galaxy CMD well confined in color. Both VI and V15 in 
the CMD lie well inside the boundaries of the theoreti- 
cal IS for Cepheids at the metallicity and reddening of 
IZwl8, while the faintest of the three variables, star V6, 
lies close to the blue edge of the predicted IS, but still 
within the Icr confidence level. We note that the V am- 
plitude of V6 seems to be slightly small compared to the 
/ amplitude (see columns 8, 10 and 11 of Table 2). The 
blue color and the slightly small V amplitude could be 
explained if the star were contaminated by an unresolved 
blue companion or by gas. ^-nd sharpness parameters 
of V6 do not seem to differ from those of a regular single 
stellar source, but the star is clearly located in a gas rich 
region (see Figure S]), thus could easily be contaminated 
by gas. The large throughput of the F606W filter to- 
wards emitting gas wavelengths along with V6 intrinsic 
faintness could thus have conspired to render the star 
V — I color bluer than normal. 

In summary, the different diagnostics all consistently 
suggest a classification as Classical Cepheids for VI, V6 
and V15. However, of the three stars only V6 (P=8.71 
days) has periodicity similar to the typical period ranges 
commonly observed for Classical Cepheids in and out- 
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side the Milky Way. This star also safely falls below the 
period of P'^ 10 days, where a break has been suggested 
to occur in the Period — Lu minosity relation o f Classical 
Cepheids (see, e.g. iN geow fc Kanbur 2004; N geow et alJ 
[200l ISa^ge et alJl2004[ ). VI and V15 instead have 
rather long periods, which, so far, have been observed 
only very r arely among Clas sical Cepheids. In a re- 
cent paper, iBird et all ()2009[ ) list the very long period 
Classical Cepheids known so f ar, which th e y nam e ultra- 
long period (ULP) Cepheids. IBird et all (2009) list in- 
cludes a total number of 18 ULP Cepheids with peri- 
ods in the range from 83.0 to 210.4 days, distributed 
over five galaxies, namely, the Large and Small Magel- 
lanic Clouds, NGC300, NGC6822, and IZwl8, referring 
to preliminary results for VI and V15 that we published 
(jFiorentino et al. 2008) . Inspection o f the literature light 
curves compiled by iBird et al.l ([2009') reveals that those 
of the ULP Cepheids seem to be more sinusoidal and 
to have a stronger dependence on wavelength than their 
classical "short" period counterparts. In most of the 
cases the / band light curves show a sort of plateau rather 
than the peak at maximum light observed in the U-band 
light curves. However, the amplitude ratios, Aj/Ay, do 
not seem to differ significantly from those of Classical 
Cepheids in standard period ranges (see last column of 
Table 2). 

As already noted in IBird et all (|2009f) , the "theoreti- 
cal" interpretation of ULP Cepheids like the bright and 
long period Cepheids detected in IZwl8 poses some prob- 
lems. In fact the current set of nonlinear convective pul- 
sation models fails to predict pulsation periods longer 
than about 90 days at this metal content for masses in 
the typical range of Classical Cepheids (3 <M< 13 M0; 
see discussion in MIO) and current evolutionary calcula- 
tions do not predict the existence of a blue loop at so high 
luminosity levels (see discussion below). Nevertheless, we 
emphasize the importance of studying Classical Cepheids 
at the very low metallicity of IZwl8 to constrain the still 
unexplored properties of the period-luminosity (PL) re- 
lation in very poor metallicity regimes. This topic is 
addressed in MIO, where we have combined the new pul- 
sation models computed to study the IZwl8 Cepheids 
with previous sets of models at higher metal abundance 
(|Marconi et al.l l2005l . , and references therein) to con- 
strain the metallicity dependence of the Cepheid PL re- 
lation in the low metallicity regime (Z from 0.0004 to 
0.008). 

The other two confirmed variables in IZwl8, V4 and 
V7, have both rather red colors. Indeed, V7 is the red- 
dest of our confirmed variables with {V — I) =1.71 mag. 
The star lies about 3 magnitudes above the RGB Tip of 
IZwl8 in the expected region of TP-AGB stars. The po- 
sition in the CMD, the long period and the shape of the 
light curve suggest that the star could be a long period 
variable (LPV) of Mira type. 

The fifth variable, star V4, is one of the brightest and 
reddest objects in the IZwl8 CMD. Only a portion of 
the star light curve is covered by our observations, thus 
both period and amplitudes of the light variation are very 
uncertain. Although the brightness and shape of the light 
curve could suggest the classification of V4 as a Cepheid, 
the very red color places the star well outside the typical 
color range of Classical Cepheids. The uncertainty of the 
period, which anyway appears to be extraordinarily long 



for a Classical Cepheid, and the rather red color, which 
seems intrinsic and not due to contamination by a red 
companion, suggest to adopt a conservative classification 
as LPV, and to not use the star to estimate the distance 
to IZwl8. 

The theoretical evolutionary tracks shown in Figure 
[2] span from the post-Main Sequence up to the early 
Asymptotic Giant Branch evolutionary phase, thus cov- 
ering both the phase of central Helium and shell Hy- 
drogen burnings. In this phase an intermediate-mass 
star crosses the pulsation IS while performing the blue- 
ward excursion usually called "blue loop" , and thus can 
be observed as Classical Cepheid. According to the 
Padua94 models, the evolutionary track that within the 
uncertainties most closely reproduces the range in mag- 
nitude and color spanned by the Cepheid with P=8.71 
days (26.2 < / < 26.8 mag, and 0.5 < ^ - / < 
0.7 mag, respectively) corresponds to M = 6Mq (red 
line in Figure [5]). The Padua tracks include a mild 
overshooting resulting in blue loops brighter by about 
0.25 dex than tracks withouth overshooting such as the 
ones predicted by canonical evolutionary scenarios [see 
also http://astro'.df.unipi.it/SAA/PEL/ZO.html. In the 
canonical assumption, the same luminosity level is at- 
tained for higher stellar masses closer to M = 6.5M0. 

On the other hand, the brighter Cepheids should have 
masses larger than 15 Mq. However, the 20 M0 track 
crosses the IS in the very fast phase of central Helium ex- 
haustion, and the coarse distribution in mass of models 
between 15 and 20 Mq does not allow us to distinguish 
the largest mass that performs the "blue loop" within the 
IS. The long period Cepheids we have detected in IZwl8 
seem to indicate that massive stars at the metallicity of 
IZwl8 (Z=0.0004) may cross the instability strip long 
enough to be observed. Whether this occurs during the 
"blue loop" or not should be further investigated from 
a theoretical point of view by both stellar evolution and 
pulsation models. Blue loops are not predicted for mas- 
sive stars by cur rent stellar models, but are invoked by 
IBird et al.l (|2009[ ) for the case of HV829, proposed to be a 
"second crossing" ULP Cepheid in the SMC (Z=0.004). 

5. THE DISTANCE TO IZW18 FROM CLASSICAL 
CEPHEIDS 

In the following we apply theoretical Wesenheit (WPL) 
relations, based on the new pulsation models computed 
by MIO, to the Classical Cepheid with P = 8.71 d (V6), 
to constrain the distance modulus of IZwl8. We do not 
apply the same relation to the Cepheids with longer pe- 
riods because pulsation models fail to predict such very 
long periods at this chemical composition. However, ex- 
trapolation of the predicted WPL relations to the very 
long periods of VI and V15 seems to well reproduce the 
data (see MIO for details). Moreover, given the small 
number of Cepheids identified in IZwl8, it is not safe to 
use the PL relations, which are known to provide reli- 
able results only when applied to statistically significant 
samples, and we prefer to use instead pulsation relations 
that can be applied to individual Cepheids, such as the 
PLC and the WPL (Madorc 1982) relations. In particu- 
lar, the latter has the ad vantage of being reddening -free 
by definition (see lMadorelll982l : iCaputo et al.ll2000i . . for 
details). 

A further method to measure the distance modulus 
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of IZwl8, based on our new pulsation mo dels, is the 
light curve model fitting technique fsee e.g..lBonq et alj 
200l 'Mar coni fc Clementinil [20051: IKeller fc Woodll200a 
Natalc ct alJ l2008L and reference therein). In the follow- 
ing subsections we will discuss and use these two meth- 
ods, separately, and will then co mpare their results with 
the distance to IZwl8 derived bv lAloisi et al.l ()2007[ ) from 
the RGB Tip method. 

5.1. Theoretical Wesenheit relations 

Updated reddening free theoretical (/, V — I) WPL re- 
lations have been derived from our new theoretical pul- 
sation models with Z=0.0004 and Y=0.24 (MIO) using 
both a canonical and a non-canonical (overluminous by 
about 0.25 dex in logi/L©) assumption on the evolu- 
tionary Mass-Luminosity relation (see MIO for details): 

< ' Mi > -IM{V 

/)=-2.58(±0.12)-3.49(±0.03)logP-h0.84(±0.1Q)logL/Lc 
with a— 0.12 mag and log L/Lc'=0-25 dex, if we use over- 
louminous models, thus including both canonical and 
overluminous solutions. 

By applying this new WPL relation to V6 we ob- 
tain /xo=31.4±Q.3 mag (D=19.0l:2:5 ^P^) in the canon- 
ical assumption, and fj,o—31.2 ±0.3 mag (D=17.4l2'2 
Mpc) from overluminous models. The distance mod- 
uli obtained by the WPL relations bracket the RBG 
Tip modulus of IZwlS obtained by lAloisi et all ()2007f ) 
(/xo=31.30±0.17 mag), and all three values are consis- 
tent within the errors. 

5.2. Light curve model fitting 

The nonlinearity of our pulsation code allows us to 
predict the variation of all the relevant physical quan- 
tities of a variable star along the pulsation cycle and, 
in particular, to model the shape and morphology of 
the light curves of individual Cepheids. This capabil- 
ity can be used to theoretically fit the observed light 
curves, thus providing another estimate of the stellar pa- 
rameters and the corresponding distance modulus of the 
host galaxy. This technique was already applied with 
success to both RR Lyrae stars and Classical Cep heids 
(see e.g-.lBono ct al. 2 0"ol iMarconi fc Clementinil 120051: 
IKeller &: WoodI i2006i: iNatale et al.l 120081 and"Teference 
therein). Here, we have applied the method to V6, using 
the new pulsation models by MIO. In particular, we ex- 
plored the mass range from 5.9 to 6.8 M©, as suggested 
by the comparison with different sets of canonical and 
non-canonical evolutionary models, and computed isope- 
riodic sequences of pulsation models by varying both the 
luminosity and the effective temperature. 

The results of the fitting procedure are shown in Figure 
[6l where we compare the observed light curves (filled cir- 
cles) to our best fit model (solid lines) , that is the model 
minimizing the mean residuals with respect to the ob- 
served curves in both V and / bands. The best fit model 
has input parameters: M=6.5 Mq, logi/L© =3.68 dex 
and Tg =5855 K, corresponding to apparent distance 
moduli of 31.5 ± 0.2 mag in both bands, where the un- 
certainty includes both the internal dispersion of the fit 
and the observational errors. Therefore, by correcting 
the obtained /zy and /z/ values for the assumed Galactic 
foreground extinction we find fio=31A ± 0.2 mag. The 
distance modulus derived with the model fitting tech- 



nique is in good agreement with the results obtained from 
the application of the WPL relations and the RGB Tip 
method. 

Taking into account all the pulsation models that fit 
the V, I light curves within a mean residual < 0.13 mag, 
we can estimate the theoretical uncertainties on the in- 
trinsic stellar parameters of the best fit model. We 
find that luminosities from 3.64 to 3.69 dex in log L/Lq 
and effective temperatures from 5800 to 6000 K are al- 
lowed, for our selection of masses (5.9 to 6.8 Mq). These 
masses, luminosities, and effective temperatures are in 
good agreement, within the uncertainties, with the val- 
ues predicted by stellar evolution models (see discussion 
in Section 4). In particular, our result that the best fit 
model corresponds to a stellar mass of 6.5 Mq seems to 
support a canonical evolutionary scenario. 

6. DISCUSSION AND CONCLUSIONS 

We have identified and derived periods for five variable 
stars in the blue compact dwarf galaxy IZwl8, among a 
list of 34 candidates, based on time series images of the 
galaxy obtained with the Advanced Camera for Surveys 
on board the HST. Three of the variable stars are Clas- 
sical Cepheids, the remaining 2 are very red objects that 
we tentatively classify as long period variables. 

The distance to IZwl8 was estimated from the Clas- 
sical Cepheids using theoretical V, I Wesenheit relations 
specifically computed for the very low metal abundance 
of the galaxy, and by theoretically fitting the observed 
light curves of the 8.71 days Cepheid. The Wesenheit 
relations provide an average distance modulus for IZwlS 
of 31.4 ±0.3 mag (0=19.0^2:5 ^pc) adopting a canoni- 
cal ML relation and of 31.2 ±0.3 mag (D=17.4l|^ Mpc) 
when a overluminosity of 0.25 dex (in log I//Lc) is as- 
sumed for each given mass. The theoretical modeling of 
the 8.71 days Cepheid light curves provides: 31.4 ±0.2 
mag (D=19.0ti-j Mpc). 

The Cepheid-based distance estimates are fully consis- 
tent within the errors wit h the distanc e modulus of 31.30 
± 0.17 mag, obtained bv lAloisi et al.l (1200 7) from a to- 
tally independent distance indicator, namely, the RGB 
Tip of IZwl8. 

Classical Cepheids and the RGB Tip of IZwlS both 
consistently indicate that the galaxy is much further 
away than previously thought, in turn supporting the 
suggestion by Aloisi et al. (2007) that IZwlS is at least 
as old as about 2 Gyrs. 

The number of Classical Cepheids we were able to iden- 
tify in IZwlS is undoubtedly very small, compared to the 
galaxy richness in gas and blue stars. Also noticeable is 
the lack in IZwlS of Classical Cepheids with periods in- 
termediate between the 10 and 100 days. As a matter 
of fact the CMD of IZwlS does not appear to containn 
many stars within and/or close to the theoretical bound- 
aries of the IS (see Figure [5]). As an additional check we 
searched for variability all the stars that lie in this region 
of the CMD, but did not succeed to identify any further 
reliable candidate. 

The gap in periods, as well as in luminosities, of the 
Cepheids in IZwl S could be related, as discussed in 
lAloisi et"al] (|2007D . to the starbursting nature of IZwlS. 
In fact, although very active in very recent epochs, 
the galaxy may lack a significant Star Formation (SF) 
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TABLE 1 
Log of the observations 



Dates 


Toxp 


HJD'' 


Filter 


Program ID 


UT 


sec 


(-2450000) 







2003- 


05 


26 


8100 


2786.129397 


F555W GO 9400 


2003- 


05 


29 


8100 


2788.995991 


F814W " 


2003- 


05 


30 


11120 


2790.301071 


F555W " 


2003- 


05 


31 


8100 


2790.995864 


F814W " 


2003- 


06 


01 


8100 


2791.995813 


F555W " 


2003- 


06 


03 


8100 


2793.795690 


F555W " 


2003- 


06 


05 


8100 


2795.928907 


F814W " 


2003- 


06 


06 


8100 


2796.879372 


F555W 


2005- 


10 


02 


2140/2140 


3646.074922/3646.200173 


F606W/F814W GO 10586 


2005- 


10 


29 


2140 


3673.354491 


F606W 


2005- 


11 


04 


2140/2196 


3679.215186/3679.282332 


F606W/F814W 


2005- 


11 


08 


2140/2196 


3683.344260/3683.411429 


F606W/F814W 


2005- 


11 


10 


2140/2196 


3685.208977/3685.276123 


F606W/F814W 


2005- 


11 


13 


2140/2196 


3688.11391273688.189299 


F606W/F814W 


2005- 


11 


18 


2140/2196 


3692.577572/3692.653076 


F606W/F814W 


2005- 


12 


03 


2140/2196 


3708.326761/3708.381510 


F606W/F814W 


2005- 


12 


10 


2140/2196 


3714.535979/3714.649539 


F606W/F814W 


2005- 


12 


20 


2140/2196 


3725.390374/3725.456858 


F606W/F814W 


2005- 


12 


25 


2140/2196 


3730.251196/3730.318190 


F606W/F814W 


2006- 


01 


01 


2140/2196 


3737.110567/3737.177722 


F606W/F814W 


2006- 


01 


06 


2008/2056 


3742.035879/3742. 103508 


F606W/F814W 



" HJD is the Helicocentric Julian Day of observation at mid exposure. 



at the epochs when the missing Cepheids should have 
for med. Indeed , once we scale the SF history derived 
by lAloisi et al.l (|1999[ ) from their assumed distance of 
10 Mpc for IZwlS to the 18.2 Mpc distance derived by 
lAloisi et al.l ()2007D , we expect few stars in the mass range 
from 6 to 20 Mq (corresponding to pulsation periods 
from 10 to 120 days) currently into the IS. A detailed 
derivation of the SF history of IZwlS from the new data 
is in progress (Annibali et al. 2010, in preparation). 
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TABLE 2 

Identification and properties of the confirmed variable stars in IZw18. 









P 


Epoch(max) 














Name'' 


(2000.0) 


(2000.0) 


(days) 


(-2450000) 


Type"^ 


{V} 


Ay 


{1} 


A? 


Ai/A^ 


VI 


9:34:1.94 


55:14:25.87 


130.3 


3726.00 


CC 


23.96 


0.91 


23.00 


0.65 


0.71 


V4 


9:34:2.30 


55:14:24.17 


196.0 


2777.00 


LPV 


23.41 


>0.29 


21.87 


>0.15 


0.52 


V6 


9:34:2.07 


55:14:21.23 


8.71 


3685.15 


CC 


27.11 


0.85 


26.48 


0.63 


0.74 


V7 


9:34:2.33 


55:14:20.49 


106.0 


3662.00 


LPV 


25.96 


>0.49 


24.25 


>0.25 


0.54 


V15 


9:34:0.46 


55:14:31.31 


125.0 


3715.055 


CC 


23.65 


0.54 


22.68 


0.41 


0.76 



" Variable stars were given increasing numbers with increasing the distance from the galaxy center which was set at 
a = 9''34'"02M0 and 6 = 55°14'26".81 (J2000.0). 

Coordinates have been measured on the HST image j9ej01010_drz.fits. 
° CC = Classical Cepheid; LPV — long period variable. 

Ay and A/ indicate the amplitude of the light variation in the V and / bands, respectively. 



TABLE 3 

Identification and properties of candidate variable stars in IZw18 



Name* 


a 

(2000.0) 


S 

(2000.0) 


{V) 


AF606W 




(I) 




AF814W 


V2 


9:34:01.76 


55; 


;14:24.75 


28.17±0.10 


0.90 


28, 


,16±0. 


,16 


0.86 


V3 


9:34:02.19 


55; 


;14;24.03 


26.43±0.04 


0.71 


26, 


.35±0, 


10 


0.64 


V5 


9:34:02.31 


55; 


;14;24.25 


26.15±0.04 


0.65 


25, 


.75±0, 


16 


0.70 


V8 


9:34:01.35 


55; 


;14;22.08 


26.50±0.03 


0.27 


24, 


.85±0, 


,03 


0.30 


V9 


9:34:02.65 


55; 


;14;21.37 


26.79±0.11 


0.90 


25, 


.83±0, 


16 


0.81 


VIO 


9:34:02.59 


55; 


;14;20.78 


26.39±0.03 


0.50 


26, 


.49±0, 


,07 


0.51 


Vll 


9:34:02.67 


55; 


;14;35.87 


28.38±0.06 


1.46 


27, 


.80±0, 


11 


0.95 


V12 


9:34:03.24 


55; 


;14;29.00 


28.79±0.08 


0.62 


28, 


.31±0, 


17 


0.35 


VIS 


9:34:03.06 


55; 


;14:20.61 


27.73±0.04 


0.66 


26, 


.59±0. 


,07 


0.71 


V14 


9:34:03.08 


55; 


;14:20.13 


26.93±0.03 


0.25 


25, 


.07±0. 


,05 


0.66 


V16 


9:34:02.74 


55; 


;14:42.78 


27.17±0.03 


0.87 


25, 


.18±0. 


,05 


1.24 


V17 


9:34:04.05 


55; 


;14:32.45 


28.13±0.05 


0.81 


26, 


.77±0. 


,07 


0.81 


VIS 


9:33:59.67 


55; 


;14:21.82 


27.41±0.06 


1.11 


26, 


.77±0. 


,13 


1.02 


V19 


9:34:00.53 


55; 


;14:43.74 


27.94±0.04 


0.98 


28, 


,22±0. 


10 


0.88 


V20 


9:34:00.19 


55; 


;14;42.60 


27.36±0.03 


0.69 


26, 


.89±0, 


,05 


0.67 


V21 


9:34:00.38 


55; 


:14:44.26 


27.60±0.03 


1.02 


27, 


.70±0, 


11 


1.45 


V22 


9:34:00.14 


55; 


:14:43.05 


26.21±0.03 


0.34 


26, 


.13±0, 


,07 


0.66 


V23 


9:33:59.78 


55; 


;14;43.99 


24.98±0.03 


0.40 


23, 


.62±0, 


,04 


0.20 


V24 


9:34:02.25 


55; 


;14;52.83 


27.50±0.04 


0.35 


27, 


.08±0, 


13 


0.35 


V25 


9:33:59.75 


55; 


;14;44.04 


24.16±0.05 


0.74 


23, 


.78±0, 


12 


1.00 


V26 


9:33:59.68 


55; 


;14;43.64 


28.32±0.06 


1.12 


27, 


.81±0. 


12 


0.93 


V27 


9:33:59.49 


55; 


;14:44.08 


27.54±0.05 


0.71 


26, 


.43±0. 


,05 


0.71 


V28 


9:33:58.60 


55; 


;14:32.22 


27.56±0.06 


0.90 


26, 


.70±0. 


,12 


1.15 


V29 


9:33:58.77 


55; 


;14:38.00 


27.42±0.04 


0.80 


26, 


.25±0. 


,06 


1.25 


V30 


9:33:58.51 


55; 


;14:38.07 


27.60±0.05 


0.83 


26, 


.68±0. 


,11 


0.75 


V31 


9:33:59.22 


55; 


;14:48.50 


26.69±0.05 


0.30 


25, 


.45±0. 


,05 


0.44 


V32 


9:33:59.46 


55; 


;14:51.53 


28.11±0.06 


1.12 


26, 


.78±0. 


,09 


1.26 


V33 


9:33:58.65 


55; 


; 14:42.92 


27.54±0.04 


0.62 


26, 


.21±0. 


,05 


1.09 


V34 


9:33:58.87 


55; 


;14:49.45 


27.84±0.05 


0.70 


26, 


.97±0. 


,10 


1.69 



" Variables were given increasing numbers with increasing the distance from the galaxy center. 
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TABLE 4 

V, / PHOTOMETRY OF THE CONFIRMED VARIABLE STARS IN IZwl8 



Star V6 - Classical Cepheid - P=8.71 days 



HJD 
(-2450000) 




V 




HJD 
(-2450000) 




I 




2786.129397 


27. 


33±0, 


,08 


2788.995991 


26, 


,29±0, 


,08 


2790.301071 


26. 


,95±0, 


,06 


2790.995864 


26, 


,53±0, 


,07 


2791.995813 


27. 


,44±0, 


,06 


2795.928907 


26, 


,05±0, 


,09 


2793.795690 


27. 


,51±0, 


,10 


3646.200173 


26, 


,63±0, 


,06 


2796.879372 


26, 


,64±0, 


,04 


3679.282332 


26, 


,64±0, 


,07 


3646.074922 


27, 


,37±0, 


,07 


3683.411429 


26, 


,83±0, 


,06 


3673.354491 


27, 


53±0, 


,07 


3685.276123 


26, 


,29±0, 


,12 


3679.215186 


27, 


,32±0, 


,06 


3688.189299 


26, 


,48±0, 


,07 


3683.344260 


27, 


,26±0, 


,06 


3692.653076 


26, 


,33±0, 


,11 


3685.208977 


26, 


,84±0, 


,04 


3708.381510 


26, 


,85±0, 


,06 


3688.113912 


27, 


,32±0, 


,10 


3714.649539 


26, 


,58±0, 


,09 


3692.577572 


26. 


,98±0, 


,05 


3725.456858 


26, 


,63±0, 


,08 


3708.326761 


27. 


,55±0, 


.10 


3730.318190 


26, 


,22±0, 


,06 


3714.535979 


27. 


,30±0, 


.06 


3737.177722 


26, 


,28±0, 


,06 


3725.390374 


27. 


,59±0, 


.09 


3742.103508 


26, 


,84±0, 


,06 


3730.251196 


26. 


,96±0, 


.04 










3737.110567 


26. 


,63±0, 


.03 










3742.035879 


27. 


,30±0, 


.05 











A portion of Table 4 is shown here for guidance regarding its form and content. The entire catalog is available in the electronic edition of 
the journal. 
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Fig. 1. — Comparison between V magnitudes obtained from the calibration to the Johnson-Cousins system of the F606W and F555W 
magnitudes, separately. A zero point shift exists between the two V-calibrated photometries for magnitudes brighter than l/=26.5 mag 
(arrow in the figure), with the F606W V-calibrated magnitudes being on average 0.04 mag brighter tfiat the F555W ones. 



Variable stars in IZwl8. 
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V-I 

Fig. 2. — I,V — I CMD of IZwlS obtained in the present study. Points have not been corrected for Galactic extintion. FiUed red triangles 
show the five confirmed variable stars detected in IZwlS, which are plotted according to their intensity-averaged mean magnitudes and 
colors. Filled orange circles show 29 candi date variables for which we could not derive reliable periods. The solid green lines are the 
evolutionary tracks bv IFagotto et al.l {1991) for masses of 4, 5, 7, 9, 12, 15, 20, 30 and 60 Mq and Z=0.0004. The evolutionary track 
that best fits the variable at (/)= 26.48 mag and {V — I)= 0.63 mag (star V6) is shown in red, and corresponds to a 6 M©. Solid and 
dashed lines represent the theoretical boundaries of the Instability Strip (IS) computed for mixing length values of 1.5 and 2, respectively 
(see Section 4 for further details). Both evolutionary tracks and boundaries of the theoretical IS were reddened for a Galactic foreground 
extinction of E{B — y)=0.032 mag (corresponding to A^=0.106 and A/=0.062 mag , respectively). The horizontal red line at I = 27.33 
mag shows the location of IZwlS RGB Tip (T-RGB), according to lAloisi et al.l II2007I ) reddened value. 
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Fig. 3. — Image of IZwlS showing the location of all 34 candidate variables (open circles), including the five confirmed variable stars. 




Fig. 4. — Image of the main body of IZwlS showing only the location of the five confirmed variable stars (open circles). 
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Fig. 5. — V (left) and / (right) light curves of the confirmed variable stars in IZwlS. Stars are labeled according to Table 2. The star's 
periods are indicated. Open and filled circles correspond to our 2007 proprietary data, and to the archival data, respectively. 
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Fig. 6. — Theoretical model fitting of the V, I light curves of the P=8.7f days Cepheid (star V6). The stellar parameters of the theoretical 
models (solid lines) best reproducing the observed light curves are labelled. Symbols are as in Figure 5. 



